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^ : ABSTRACT 

^ ' We discuss the large radiative corrections to the masses of supersymmetric Higgs bosons 

OO . in the MSSM and in its simplest extension, the NMSSM, with particular attention paid to the 

' bounds on the lightest CP-even Higgs mass found in both models. In the case of the MSSM, these 

. corrections are found to be primarily associated with the effects top quark and stop squark loops, 

' while for extended models they also include significant contributions from Higgs and Higgsino loops. 

m ' 

^ . 1. Introduction - SUSY Higgs Bosons at Tree Level 



1.1. MSSM 



p • The Minimal Supersymmetric Standard Model[|T| (MSSM) is the simplest 

1^ . possible supersymmetric theory which symmetry breaking can reduce to the Stan- 

dard Model. It contains not only the usual fermionic sector with SU(3) xSU(2) xU(l) 
gauge bosons plus their superpartners, but also two supermultiplets of Higgs dou- 
^ ■ blets. This is the minimum Higgs sector which allows for the correct electro-weak 

breaking structure and mass couplings, because supersymmetry both forces the in- 
clusion of higgsinos, which cause anomalies which cannot be cancelled with only one 
doublet, and restricts the possible Yukawa couplings between Higgses and fermions. 

We shall use the convention that the two doublets are labelled Hi and H2, 
with H2 coupling to the up type quarks. The superpotential is then given by 

W = huQH2U + hdQHid + h.LHie + i^tHiH2 (1) 

where we have only shown one generation. When supersymmetry is broken, the 
resulting scalar potential acquires soft breaking terms of the form 

Vsoft ^^mf\(l)}\+huAuQH2U + hdAdQHid + h^A^LHie + iiBHiH2 (2) 

i 

where the sum is over all spin particles, and we use the usual convention of using 
the same label for the scalars as for the full supermultiplet. 

*talk presented at SUSY-93, Boston, March 1993 



1 



With two doublets we shall have more physical Higgs particles - two CP-even 
(labelled h for the lighter and H for the heavier), one CP-odd (or axial, labelled A), 
and one charged (labelled - in additional to the usual Goldstone modes. 

It is simple to derive the following Higgs potential from the lagrangian : 

Vmggs ^ Y^HlH.f + ^{HlH2f + {\^ + \,){H\H,){HlH2)-\i\HlH,\'' 

+ mliHlHi) + ml{HlH2) ~ml{H^ia2H2 + h.c.) (3) 

where we may use the requirements of supersymmetry to obtain 

M = M = — — A3 = — - — ^ Y ^ ' 

and the nii are functions of the unknown soft parameters of the theory. Thus we may 
calculate all of the Higgs boson masses in terms of the (known) gauge couplings and 
three unknown masses. Defining < H° >= vi, < >= t^2, tan/3 = ^, and v"^ = vl + 
we may use the relation = \{g\ + 52)^^ to remove one of the unknown masses. 
It is then conventional to reparametrise in terms of tan/? and m^, the mass of the 
CP-odd particle. Once we have done this, we find 

2 2,2 

vdn — m A ^ m,, 



''A "1" '"'W 



2 _ ^ /™2 I 2 



± Wmj^ + TO^ — 2m^TO| cos4/3 ) (5) 



where mc is the mass of the 

From this analysis, we should note primarily that the mass is larger than 
the W mass; while the h mass is less than both and mz. This forms an extremely 
firm prediction about the Higgs sector. We should note, however, that as yet we 
have not discussed radiative corrections, and are working purely with a low energy 
effective theory including only Higgs bosons. 

1.2. NMSSM 

The Next-to-Minimal Supersymmetric Model (NMSSM) was discussed in 
detail in the talk by Terry Elliott at this conference |p. This model has a Higgs 
structure identical to that of the MSSM, but with an extra singlet, and an extra 
term in the superpotential \NH1H2 + |A^^ to replace the parameter of the MSSM. 

As was also described earlier at this conference ^, the consideration of 
the perturbative bound on the coupling A gives a bound on the mass of the lightest 
Higgs boson which is rather larger than that found for the MSSMH. 



2. Radiative Corrections 

So far we have calculated masses using a low energy effective theory from 
which all heavy particles (including all superpartners) have been removed. The first 
problem with this analysis is that we have used relations defining the couplings Aj 
which are only strictly valid at some high energy scale which we call nisusy, above 
which supersymmetry is a good approximation. Below this scale, the couplings will 
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run from their starting values. This suggests a renormahsation group analysis to ac- 
count for the difference between the values of the parameters at the supersymmetry 
and Higgs (weak) scales. 

Furthermore, in decoupling the superpartners at high energy, we have ne- 
glected any threshold effects, which may come both from the fact that the su- 
perpartners are not all degenerate with mass rususy but have some spectrum, and 
from the finite one-loop diagrams involving these particles which may give non-zero 
contributions to the Higgs masses. These effects can be best evaluated by either 
simply calculating the diagrams at some high energy, or by carrying out a one-loop 
Coleman- Weinberg effective potential calculation. 

Lastly, we should also be careful that loops of light Higgs particles do not 
themselves contribute substantial radiative corrections to the mass matrices. 

We shall consider only the top and stop (and, for the charged mass, the 
sbottom) as coupling strongly to the Higgs sector, as all of the other particles 
couple through Yukawa couplings which are relatively very small (unless tan/3 >> 1 
in which case the bottom quark Yukawa is also large; we shall not consider this 
case). Further corrections may also be given by loops of Higgses and Higgsinos. We 
shall neglect the consideration of wave-function renormahsation and of gauge boson 
loops as relatively small. The former is negligible so long as the masses which we 
calculate are small relative to those of the particles in the contributing loops. 

3. Results - MSSM 

The situation in the MSSM is relatively simple, because as we have seen the 
couplings involving Higgses and Higgsinos are small (of order gl at worst), and so 
the most important effects are from top and stop loops. If we begin by considering 
only top effects, we recall that the top mass is given by mt = htvsinfi, where ht is the 
top quark Yukawa coupling. The simplest way of finding the top quark contribution 
is to use the renormahsation group (RG) equations below the supersymmetry scale 
p. These are rather long, but the largest correction^ is given by 



where t is the log of the renormahsation point. It is then trivial to solve this reduced 
form to give the result 



which has also been obtained using effective potential techniques [§. In fact, solving 
the full RG equations shows that this estimate is rather conservative so that the 
bound should be rather lower, because the effects proportional to hlgl are of the 
opposite sign. In any case the most important correction here being proportional 
to hi (and thus to mf), it becomes large very rapidly with increasing top mass. 

Typical sizes of the increase in the bound found numerically are 2, 11, 22, 
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35 GeV for a top mass of 120, 140, 160, 180 GeV and m^usy =lTeV, and for the 
simplified analytical approximation given above are around 5 to 8 GeV more. 

Next we consider the effects of squarks. These are rather messy, and to per- 
form this calculation we may either do an effective potential calculation or evaluate 
all the appropriate graphs. The results have been given in full by Ellis et a/[Q, and 
we shall not present them here, but merely mention the interesting features. Firstly, 
the terms which are not logarithms involving the renormalisation scale are all pro- 
portional to At + fj, cot i3 to some power. This is the off-diagonal term in the stop mass 
matrix (to within factors) and so these corrections are zero unless there is mixing 
between left and right-handed stops. Next, the corrections are all proportional to 
hf. Finally, as we would expect by supersymmetry, the logarithmic terms from the 
divergent diagrams exactly cancel the top corrections if the top and both stops are 
degenerate. After some algebra, one can derive the following analytic formula for 
the bound by simply maximising the formulae with respect to all the parameters: 



where we use m^^ to indicate the heavier of the stop mass eigenstates. Note that 
the nisusy dependence has been cancelled; for puposes of calculating top and stop 
loops, nisusy is effectively the stop mass. 

The finite terms terms are relatively insensitive to the stop masses except 
indirectly through the constraint that colour must not be broken, and so the stop 
mass eignvalues must be positive, which gives 2htV2iAt + ^J■cotp) < (to?^ -"t-^J; however 
they are sensitive to At and At + ^j, cot /3, and so for most arbitrarily chosen values of 
the parameters the lighter Higgs mass will not approach this bound, and in fact for 
very large At the corrections to the bound may even be negative. 

In summary, Eq. (P) is a pessimistic view of the final bound after radiative 
corrections; one would not expect it to be realised except in the rather unlikely case 
that all of the parameters are such as to make the mass as large as possible, but 
if this bound is covered by the searches at LEP then we can safely rule out the 



In addition to these corrections to the bound, large corrections also apply 
to the mass spectrum, particularly in the region where tan^ is small (the bound is 
saturated as tan/3 oo; we shall not discuss these here in detail, except to note that 
they are particularly affected by contributions from gauge loops. This is because 
the top quark decouples at a scale to* (and similarly for stop squarks), while the 
Higgs sector continues to contribute to the RG equations down to lower energies. 

4. Results - NMSSM 

Unfortunately, this rather simple analysis is drastically complicated in the 
case of the NMSSM. Apart from the obvious point that 3x3 matrices are not 
susceptible to simple diagonalisation, the spectrum is much more complicated, and 
so our RG analysis approach, which implicitly assumed that all Higgs bosons had at 
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least approximately similar masses, is invalidated. Apart from this, there are now 
large (order 1) couplings in the Higgs and Higgsino sectors. This will mean that we 
must carry out an analysis of the effects of loops of these particles. 

The effects considered in the MSSM are still of course present[0, and indeed 
the formulae for the radiative corrections due to top loops is still valid. The stop 
effects for this model have been derived elsewhere and their effects on the 
bound have been studied ITTI] , with the result that for the region of high top mass 



their effects on the bound may be substantial, of order 20GeV; however, this region 
is not the one where the bound is maximised, and at low top mass, where the bound 
reaches its highest values, these corrections are negligible. The main effect is thus 
to lift up the lowest part of the nih - curve, diluting the strong mt dependence of 
the bound. 

Now we must also consider effects of Higgs and Higgsino loops. In the event 
that we make the assumption that we have a region of parameter space in which 
all of the Higgs states are light (that is there are none which differ from our renor- 
malisation scale by more than, say, a factor of 2 or so), then we may use our RG 
approach to find the radiative corrections, and thus include all divergent diagrams. 
This leads to the predicted result that the corrections proportional to A are in fact 



quite large |]T2[; the effect on the bound is to raise it by approximately 7GeV (but 
by rather less for very large and small top mass), and this result is that previously 
presented by Terry Elliott |0. 

In the event that the Higgs spectrum includes heavy mass eigenstates, or if 
we wish to include Higgsino effects or effects from finite Higgs loops, the situation is 
much more complicated. The only technique which we can use to deal with this case 
is that of Coleman- Weinberg analysis, writing out all the tree level mass matrices 
and then using them to calculate the effective potential. Unfortunately, while for 
scalars this is quite simple, since the mass Higgs and Higgsino mass matrices reduce 
to the usual 3x3 and 2x2 blocks, there is no such simple result when we consider the 
charged and CP-odd particles, and we must consider the full 10x10 mass matrix. 
This is impossible to do analytically, and so we must carry out the differentiations 
numerically. This work is still in progress, but it seems that the radiative corrections 
are typically of order a few GeV and increase the calculated scalar masses for most 
values of the parameters, reaching larger values if r >> 1, where they may also give 
corrections of opposite sign. 

In conclusion, it is relatively simple to find the radiative corrections to the 
bound from the effects of top and stop loops, the latter being more difficult to do 
with precision; but a full analysis of the Higgs and Higgsino corrections to the bound 
awaits the calculation of the analytic effective potential corrections. Furthermore, 
this depends on arbitrary parameters (such as masses of Higgs bosons which are 
not directly involved in the bound) to such an extent that it is difficult to find 
results without making assumptions about which region of parameter space is most 
interesting. However, for a given set of parameters it is simple to find the full 
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spectrum including radiative corrections. 

5. Conclusions 

In conclusion, the effects of radiative corrections on the Higgs mass spec- 
trum are extremely large in both models considered here. In the MSSM they are 
dominated by the well-known contributions of top and stop loops, while for the 
NMSSM they also include comparably large and (usually) positive corrections from 
loops of Higgses and Higgsinos; however in the latter model it is much harder to 
make concrete predictions with regard to the bound because in general these effects 
are quite strongly dependent on the spectrum and thus on which region of the very 
large parameter space is used for the calculations. 

One question which we might ask is how dependent the results are on the 
specific model considered. Clearly, going from the minimal to next-to-minimal 
models has drastically increased the complexity and size of the radiative corrections, 
and it seems likely that introducing further singlets, triplets, or other structure 
would make things even worse. Given the size of the corrections (and difficulty in 
doing calculations because of the size of the parameter space) in the NMSSM, we 
should be rather wary of bounds in such models. 
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